Density preturbations in a finite scale factor singularity universe 
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We discuss evolution of density perturbations in cosmological models which admit finite scale 
factor singularities. After solving the matter perturbations equations we find that there exists a set 
of the parameters which admit a finite scale factor singularity in future and instantaneously recover 
matter density evolution history which are indistinguishable from the standard ACDM scenario. 
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I. INTRODUCTION 

One of the most problematic phenomena resulting from 
observations of high-redshift type la supernovae (SNIa) is 
recent accelerated expansion of the universe. Search for 
the explanation of this phenomena led physicists to many 
various possible cosmological scenarios based on differ- 
ent approaches like modifying physical expansion history 
or modifying the theory of gravity. Out of this effort 
there arised a couple of cosmological scenarios. Some of 
them admit new types of singularities which has already 
not been known, within the framework of the so-called 
standard or concordance cosmology. Finite scale factor 
singularities (FSF) are one of the types and were first 
found in Ref. Basically, it is assumed that the uni- 
verse is accelerating due to an unknown form of energy 
which phenomenologically behaves as the cosmological 
constant. More observational data Q made cosmologists 
think of an accelerating universe filled with phantom Q 
which violated all energy conditions: the null {q-\-p > 0), 
weak > and g -\- p > 0), strong 4- p > and 
g + 3p > 0), and dominant energy {g > 0, —g < p < g) {g 
is the energy density and p is the pressure). A phantom- 
driven dark energy leads to a big-rip singularity (BR, or 
type I according to [l|) in which the infinite values of the 
energy density and pressure (p, p oo) are acconipanied 
by the infinite value of the scale factor (a oo) 

The list of new types of singularities contains: a big- 
rip (BR) 0], a sudden future singularity (SFS) [Mil, 
a generalized sudden future siiigularity (GSFS), a finite 
scale factor singularity (FSF) '12] ja big-separation sin- 
gularity (BS) and a w-singularity [l5|. A weaker version 
of the Big- Rip such as a Little- Rip and a Pseudo-Rip has 
also been proposed recently [l^, [l^. In this paper we 
deal with a finite scale factor singularity. This is a weak 
singularity according to Tipler and a strong singularity 
according to Krolak [T3 |. 

In Ref. [l^ we found that there is a set of the pa- 
rameters which, within the la CL, fits the observational 
data BAO, SNIa and the shift parameter, and admits an 



FSF singularity. In this paper we deal with the prob- 
lem of growth of density perturbations in the scenario 
admitting such a singularity. 

The paper is organized as follows. In section [II] we 
present an FSF scenario. In section IIIII we present the 
expressions for the evolution of linear density perturba- 
tions of matter in general relativity, and rewrite them for 
the scenario admitting an FSF singularity. In section IIVI 
we give the results and discussion. 



II. A FINITE SCALE FACTOR SINGULARITY 
UNIVERSE 

In order to obtain an FSF singularity one should start 
with the simple framework of an Einstein- Friedmann cos- 
mology governed by the standard field equations (we as- 
sumed fiat universe) 
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Similarly like in the case of an SFS, which were tested 
against the observations in Refs. pTl - {23j . one is able to 
obtain an FSF singularity by taking the scale factor in 
the form 
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(II.3) 

with the appropriate choice of the constants S,ts, as,m,n. 
In contrast to an SFS, in order to have an accelerated 
expansion of the universe, 6 has to be positive {d > 0). 
For 1 < n < 2 we have an SFS. In order to have an 
FSF singularity instead of SFS, n has to be limited to 
< ?i< 1. 

As can be seen from pi.ip - (|II.3l) . for an FSF p diverges 
and we have a — >■ a^, p — > oo, and \p\ -> oo for t tg. 
In the model (jll.3[) , the evolution begins with a standard 
big-bang singularity at i = for a — 0, and finishes at a 
finite scale factor singularity at i = is, where a — Us = 
a(ts) is a constant. In terms of the rescaled time y, we 
have a(l) = a^. 
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The standard Friedmann limit (i.e. models without an 
FSF singularity) of (|II.3p is achieved when 5 — >■ 0; hence 
6 is called the "non-standardicity" parameter. Addition- 
ally, notwithstanding Rcf. 7], and in agreement with 
the field equations (|II.ip - pi.3l) . S can be both positive 
and negative leading to an acceleration or a deceleration 
of the universe, respectively. 

To our discussion it is important that the asymptotic 
behaviour of the scale factor (III. 31) close to a big-bang 
singularity at < = is given by a simple power-law 
obb = 2/™, simulating the behaviour of flat (fc — 0) 
barotropic fluid models with m = 2/[3(w + 1)] where 
w is barotropic index {p = wp). 

Recently, an FSF singularity scenario was confronted 
with baryon acoustic oscillations, distance to the last 
scattering surface, and SNIa [l^. It was shown that for a 
finite scale factor singularity there is an allowed value of 
TO = 2/3 within la CL, which corresponds to a dust-filled 
Einstein-de-Sitter universe in the past. It was also shown 
that an FSF singularity may happen within 2 x 10^ years 
in future in la confidence level, and its observational pre- 
dictions at the present moment of cosmic evolution can- 
not be distinguished from the predictions given by the 
standard quintessence scenario of future evolution in the 
Concordance Model [IIM [sil-lsl . 
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III. LINEAR DENSITY PERTURBATIONS 

In the linear regime the equations that govern the evo- 
lution of perturbations in a Friedmann universe consist- 
ing of more than one component constitute a complicated 
set of coupled differential equations [4l|. In this paper 
we consider the evolution of perturbations in a flat Fried- 
mann universe made up of a dust matter with the den- 
sity pm and a dark e nerg y with density p^e, and pressure 
Pde- It was shown in [30| that, in similar case, neglecting 
perturbations in dark energy one makes some particu- 
lar, unintended choice of gauge and in general that may 
lead to errorneous results for perturbations in the matter. 
Taking that into account we restrict our investigations to 
the cases where the proper wavelength of perturbations 
is much smaller than the Hubble radius and the sound 
velocity for the dark energy has a positive value of order 
of unity, while the barotropic index for the dark energy is 
a reasonable slowly varying function of the cosmic time. 
With these assumptions the dark matter perturbations 
effectively decouple from perturbations in the dark en- 
ergy and the evolution of the matter density contrast 6m 
can be described to a good approximation with the fol- 
lowing equations: 
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(III.l) 

(III.2) 
(III.3) 



0.15 0.51 ±0.11 
0.35 0.7 ±0.18 
0.55 0.75 ±0.18 
1.4 0.9 ±0.24 
3.0 1.46 ±0.29 



TABLE I: Observational data for the growth factor / from 
Lyman-a forests and galaxy redshift distortions taken from 
Refs EB^- 



where a dot denotes a derivative with respect to time 
H = a/a is the Hubble parameter. 

Since the amplitude of the anisotropy is determined by 
the typical amplitude of peculiar velocities and those, in 
the linear theory, correspond to growth rate of perturba- 
tions / = d\nS/dhia, we will write down the equation 
pil.ip in terms of this logarithmic growth rate as follows: 
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where da; = d In a and 
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Further, taking the scale factor pi.3p . we get the equation 
for the growth rate evolution in the form: 



(IH.6) 



where / denotes derivative with respect to y from pi.3p 
and 



a^n''' 



(in.7) 



where H = a'^/aQ and aq — a(yo)- 

Observational data from Lyman-a forests and galaxy red- 
shift distortions for the growth rate /, to which we fit our 
model, are given in table HI 



IV. RESULTS AND CONCLUSIONS 

The main goal of this paper is to find the fit to cur- 
rently available data for the growth of perturbations rate 
taken from Refs. [35l - l40l | (see tableHl for the cosmological 
model which admits an FSF singularity. We are search- 
ing for the fit, varying the model parameters which are 
TO, n, 5, 2/0, /o, and the last parameter /o is the present 
value of the growth rate /. We search for such a set of 
the parameters that satisfies, within la CL, BAO, SNIa, 
and the shift parameter data as well (see [13j). 
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FIG. 1: Upper left panel presents the predicted growth rate plotted against redshift for an FSF model (dashed line) and 
a ACDM model (solid line). There are also shown the observed data points (dots) for the growth rate. The upper right 
panel shows the distance-redshift relation for both FSF and ACDM models. The left and right bottom panels present the 
corresponding relative errors. 



We solve the equation (jIII.6[) numerically for a given 
set of the parameters using the standard Runge-Kutta 
method with an adaptative step size. Applying a stan- 
dard Levenberg-Marquardt method, we search for the 
minimum of the function which is of the form 

x(2;p) = 2^ -2 , (iv.i) 

1=1 i 

where: /obs and Oi are taken from the tableU] jth is calcu- 
lated by solving the equation piI.6P : p = (m, n, 5, yo, /o)- 
We find the following fit for one of the possible set of pa- 
rameters: 

yo = 0.55, (5 = 0.67, m = 0.49, n = 0.32, /o = 0.53, 

with — 0.99. For this set of parameters we evaluate 
the growth rate function /, again solving numerically eq. 
(IIII.6I) , cf . the upper left panel of figure [T] In this panel 
together with growth rate for FSF scenario, we see the 
growth rate for ACDM scenario and the measured values 
of the growth rate with their errorbars. In a bottom left 
panel we see the relative difference between the evalua- 
tion of the growth rate function, for an FSF scenario and 
for a ACDM. The discrepancy for both models is at most 
9% for z - 1. 



In the top right panel of the figure [T] we see the 
distance-redshift relation for an FSF scenario and a 
ACDM model. In the bottom right panel of the same 
plot we see a relative difference between distance-redshift 
relations for both models, which is biggest (14%) for the 
most distant values of z ~ 3.5. 

As in Refs. [l^ and [1^, the set of the parameters that 
we obtained was tested against several additional condi- 
tions, what assured, that some other physical conditions 
are satisfied which are listed bellow: 

• we assumed that, the scale factor and its first 
derivative for all times is always positive, i.e. 
a(y) > 0, and a{y) > 0; 

• a current expansion of the universe should be ac- 
celerated, i.e. d{yo) > 0; 

• time should be decaying function of z, positive red- 
shift should correspond to the past (z > for 
y < yo) and negative redshift should correspond 
to the future (z < for y > yo)- 

We conclude that for the FSF models there exists a 
set of parameters which fits the observational data for 
the growth rate and on the other hand satisfies, within 
1(7 CL, the data for BAO, SNIa and the shift parameter. 
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Thus we proved that current observations are incapable Center grant No N N202 3269 40 (years 2011-2013). 
of ruhng out FSF models of the expanding universe. 

V. ACKNOWLEDGEMENTS 

We acknowledge the support of the National Science 



[1] S. Nojiri, S.D. Odintsov and S. Tsujikawa, Phys. Rev. D [19 

71,063004 (2005). 
[2] S. Perlmutter et al., Astroph. J. 517, (1999) 565; A. G. [20 

Riess et al., Astron. J. 116, 1009 (1998); A.G. Riess et 

al., Astroph. J. 560, 49 (2001). [21 
[3] J.L. Tonry et al., Astroph,. J. 594, 1 (2003); M. Tegmark 

et al., Phys. Rev. D69, 103501 (2004); R.A. Knop et al., [22 

Astrophys. J. 598, 102 (2003). 
[4] R.R. Caldwell, Phys. Lett. B 545, 23 (2002); M.P. [23 

Dcibrowski, T. Stachowiak and M. Szydlowski, Phys. 

Rev. D 68, 103519 (2003); P.H. Frampton, Phys. Lett. [24 

B 562 (2003), 139; H. Stefancic, Phys. Lett. B586, 5 

(2004); E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. [25 

Rev. D70, 043539 (2004); S. Nojiri and S.D. Odintsov, [26 

Phys. Lett. B595, 1 (2004). [27 
[5] R.R. Caldwell, M. Kamionkowski, and N.N. Weinberg, 

Phys. Rev. Lett. 91, 071301 (2003). [28 
[6] V. Sahni and Yu.V. Shtanov, Class. Quantum Grav. 19, 

LlOl (2002). [29 
[7] J.D. Barrow, Class. Quantum Grav. 21, L79 (2004). 
[8] J.D. Barrow and Ch. Tsagas, Class. Quantum Grav. 22, [30 

1563 (2005). [31 
[9] S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 103522 

(2004). [32 
[10] K. Bamba, S. Nojiri, and S. D. Odintsov, J. Cosmol. 

Astropart. Phys. 10 (2008) 045. [33 
[11] M.P. D9,browski, Phys. Rev. D71, 103505 (2005). [34 
[12] M.P. Dabrowski and T. Denkiewicz, AIP Conference Pro- [35 

ceedings 1241, 561 (2010). [36 
[13] T. Denkiewicz, arXiv: 11 12.5447, 
[14] L. Fernandez-Jambrina and R. Lazkoz, Phys. Rev. [37 

D70, 121503(R) (2004); L. Fernandez-Jambrina and R. [38 

Lazkoz, Phys. Rev. D74, 064030 (2006). 
[15] M.P. Dabrowski and T. Denkiewicz, Phys. Rev. D79, [39 

063521 (2009). 

[16] P. H. Frampton, K. J. Ludwick, R. J. Scherrer, Phys. Rev. [40 
D84, 063003 (2011). [41 

[17] P. H. Frampton, K. J. Ludwick, R. J. Scherrer, 
[axXrv :1112. 2964^1 . 

[18] A. A. Starobinsky, Phys. Lett. 91B, 99 (1980). 



M. R. Setare and E. N. Saridakis, Phys. Lett. B671, 331 
(2009). 

V. Sahni and Yu. Shtanov, Phys. Rev. D 71, 084018 

(2005) . 

M. P. Dabrowski, T. Denkiewicz, M. A. Hendry, Phys. 
Rev. D75, 123524 (2007). 

H. Ghodsi, M. A. Hendry, M. P. Dabrowski, T. 
Denkiewicz, MNRAS, 414: 15171525 (2011). 
Tomasz Denkiewicz, Mariusz P. Dcabrowski, Hoda Gh- 
odsi, Martin A. Hendry. [arXiv:1201 . 6661 
M. Chevallier and D. Polarski, Int. J. Mod. Phys. DIO, 
213 (2001). 

E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003). 

E. V. Linder, Phys. Rev. D72, 043529 (2005). 

T. Koivisto and D. F. Mota, Phys. Rev. D73, 083502 

(2006) . 

R. Caldwell, A. Cooray and A. Melchiorri, Phys. Rev. 
D76, 023507 (2007). 

P. Zhang, M. Liguori, R. Bean and S. Dodelson, Phys. 
Rev. Lett. 99, 141302 (2007). 

A. J. Christopherson Phys. Rev. D82, 083515 (2010). 
L. Amendola, M. Kunz and D. Sapone, JCAP 0804, 013 
(2008). 

C. Di Porto and L. Amendola, Phys. Rev. D77, 083508 
(2008). 

W. Hu and L Sawicki, Phys. Rev. D76, 104043 (2007). 

E. V. Linder, Phys. Rev. D79, 063519 (2009). 

L. Guzzo et al, Nature 451, 541 (2008). 

M. CoUess et al, Mont. Not. R. Astron. Soc. 328, 1039 

(2001). 

M. Tegmark el al, Phys. Rev. D 74, 123507 (2006). 
N.P. Ross et al, Mont. Not. R. Astron. Soc. 381, 573 

(2007) . 

J. da Angela et al, Mont. Not. R. Astron. Soc. 383, 565 

(2008) . 

P. McDonald et al, Astrophys. J. 635, 761 (2005). 

V. F. Mukhanov, H. A. Feldman and R. H. Branden- 

berger, Phys. Rep. 205, 215 (1992). 



